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Abstract 



Ultraviolet renormalons, contrary to their infrared counterparts, lead to a universal con- 
tribution to the large-order behaviour of perturbative expansions in QCD. In this letter, 
we determine nature of the leading ultraviolet renormalon singularity for the inclusive 
hadroproduction cross section in e + e~ annihilation, for hadronic r decays and the mo- 
ments of deep-inelastic scattering structure functions. We comment on the relevance of 
ultraviolet renormalons to estimates in low orders of perturbation theory. 



1. Perturbative expansions of physical quantities in QCD are divergent, and as such 
they are related to a measurement only through the additional assumption that the se- 
ries is asymptotic. Of the two known mechanisms that cause divergence of the series, 
instantons Jl], Q] are unimportant in large orders and the asymptotic behaviour of the se- 
ries expansion is determined by renormalons p, ^| . According to whether the divergence 
arises from regions of small or large momentum in internal integrations, renormalons are 
classified as infrared (IR) or ultraviolet (UV). The growth of perturbative coefficients 
due to IR renormalons depends on whether one considers, for example, deep-inelastic 
structure functions, hadronic event shape observables or the hadronic total cross section 
in e + e~ annihilation and is closely related to power corrections to these observables. 
For this reason they have recently been studied intensely ||. UV renormalons are of- 
ten regarded to be more complicated and have mostly been left aside, perhaps because 
they are Borel summable. However, the structure of UV renormalons is universal in 
the sense that it depends only on the theory under consideration, in our case QCD. 
Moreover, for some quantities of interest, such as the hadronic total cross section in 
e + e~ annihilation, the sign-alternating UV renormalon behaviour determines the actual 
large-order behaviour of the series expansion. In this letter we determine the universal 
ultraviolet renormalon asymptotics in QCD and apply it to a number of observables of 
phenomenological interest. 

The universality of UV renormalons was recognised by Parisi ||, who noted that 
UV renormalons could be compensated by adding higher- dimension operators to the 
Lagrangian just as logarithmic UV divergences can be compensated by dimension-four 
counterterms. To state Parisi's hypothesis precisely, we consider a quantity R(a s ), ex- 
panded as 



The dependence on external momenta is not indicated and a s denotes the coupling 
renormalised at a scale //. The UV renormalons produce poles at t = m/ (3$ < in the 
Borel transform B[R](t) = J2n=o r nt n /n\, where m is a positive integer and (3q = —b/ (47r) 
(6 = 11 — 2Nf/3) the first coefficient of the /^-function. It follows that the integral^ 



is complex, for t > —1//3q, and its imaginary part is unambiguously related to UV 
renormalon singularities, and therefore asymptotic behaviour, and vice versa. In the fol- 
lowing, we will be concerned explicitly only with the leading UV renormalon singularity 

1 Since UV renormalons produce sign-alternating factorial behaviour, they do not lead to ambiguities 
in the usual Borel integral. Because the consideration of ambiguities will simplify the renormalisation 
group considerations below, we define the integral parallel to the negative axis. The imaginary parts 
created by UV renormalon poles are now exponentially large in a s . This need not bother us, because 
a s could be considered negative without any change in our derivation. 
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at t — 1/Po- The hypothesis of Parisi ||, for this situation, states that 

lmI[R](a s ,p k ) = \Y, e ' imas)) ^ llPlc ^s)Ro^Pk). (3) 

A* i 

Thus, the leading UV renormalon behaviour is determined by single zero-momentum 
insertions of dimension-six operators^] Oi into the Green function from which the quantity 
R is derived. The important point here is that the coefficient function Cj is universal - 
it is independent of the external momenta pk and the quantity R - which is related to 
the fact that UV renormalons arise from loop momentum regions much larger than any 
external scale. The dimension-six operators may be thought of as an additional term 

A£ = -- 2 £ c -VW.«.)) a7 A//3g Cl {a s ) O t (4) 

A* i 

in the QCD Lagrangian with coefficients such that for any R the imaginary part of 
I[R] is compensated by the additional contribution to R from AC By comparing the 
renormalisation group equations for I[R] and Ro i one derives that 

d CM = ^r4Q(a s ), (5) 



da s 2f3(a s 

where j(cx a ) is the anomalous dimension matrix pertaining to the dimension-six operators 
Oi. The solution to this equation determines completely the o^-dependence of the Cj. 
This in turn allows us to compute the Borel transform B[R](t) in the vicinity of the 
singular point t = 1/ (3q up to an over-all constant. Finally, the nature of the singularity 
determines the large-order behaviour. These manipulations can be summarised in the 
substitution rule 

e -l/(A,a.) /^ oC gA _ I V- pn n] n -X ft n+l (g) 

to obtain the leading asymptotic behaviour of R from the (^-dependence of Im I[R]. 

The leading UV renormalon divergence has been subjected to detailed diagrammatic 
study in abelian models []7|]-[]ID|], which confirmed the general structure of (^). In [0 



four-fermion operators were identified as sources of leading behaviour in QED. In [10 
(d) was shown to be valid in QED to all orders in the 1/Nf expansion and it was 
indicated how one could understand the restoration of the full non-abelian (3q starting 
from a 1/Nf expansion of the non-abelian theory. Both J7], [To|] corroborated the earlier 
expectation |11|-|[T3[ that while (|5]) determines the a s -dependence of Cj, the integration 



constants, which are related to the constants that specify the over-all normalisations 
K~i of the UV renormalon asymptotic behaviour, could not be calculated systematically, 
except in expansions like 1/Nf. Once this is realized, all exact obtainable information 
can be obtained solely by solving the renormalisation group equations for C{ above. This 

2 The fact that higher (than four) dimensional operators start at dimension six determines that the 
position of the leading UV singularity is indeed at t = 1/(3q. 
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problem has already been solved for Heavy Quark Effective Theory || and QED ]7|, |Kj . 
In the following sections we treat the slightly more complicated case of QCD. That is, 
for 

r n n =°° (3%n\ n A//3 ° ]T K t n Si (1 + 0(1/ n)) (7) 

i 

we determine the so far unknown constants 5j. Note that the 1/n corrections to the 
asymptotic behaviour are in principle calculable as well. This would require the two-loop 
anomalous dimension matrix as well as the one-loop corrections to the Green functions 
with operator insertion, Rd- The constants Ki remain unknown. However, because of 
the universality of coefficients functions Ci, the ratio of K^s for different observables is 
calculable. 



2. In this section we determine the leading UV renormalon behaviour of current 
correlation functions 

• / d 4 xe^ (0\T(j,(x)U0))\0) = (q,q v - g, u q 2 ) U(q 2 ) (8) 

in massless QCD with Nf flavours. We will consider colour- singlet vector and axial- vector 
currents and let them also be flavour-singlets. In the following section, we generalise to 
the real case where flavour symmetry is broken by electric or axial charges or the current 
is flavour non-diagonal. Thus, in expressions like (ipMifj) a sum over flavour, colour and 
spinor indices is implied and M is a matrix in colour and spinor space, but unity in 
flavour space. 

To account for the external currents, we introduce two U(l) background fields that 
couple to the vector and axial-vector current and consider the Lagrangian 

£QCD+jv v »+jA a n- (9) 

The corresponding field strength tensors are defined as usual and denoted by and 
H^, respectively, such that d^F^ = j v and d^H^ = f A . A basis of dimension-six 
operators is then given by 



o 1 = 




o 2 = 




o 3 = 




o A = 


(^75^(^75^) 


o 5 = 


~ J ABC tr^p (J 



(10) 



(11) 

y 

0, = - 2 $1^) d v F»* 7 = \ (^75<A) d v W (12) 
9 9 

°s = ^ 9^ d p F p , 9 = \ d u H^ d p H m . (13) 
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We need not consider gauge-variant operators and those that vanish by the equations 
of motion as their zero-momentum insertions into the above correlation functions do 
not contribute. The over-all factors l/g k have been inserted for convenience. Note that 
we did not include four-fermion operators of scalar, pseudo-scalar or tensor type. They 
can not be generated in massless QCD, because the number of Dirac matrices on any 
fermion line that connects to an external fermion in a four-point function is always 
odd. The coefficients Cj corresponding to these operators therefore vanish exactly. It 
is straightforward to compute the leading-order anomalous dimension matrix 7, defined 
such that the renormalised operators satisfy 



0. 



(14) 



We find 7 = 7 (1) a s /(47r) with 3 



7 



(i) 



(A B 

755 

\o c 



(15) 



and (C F = (Nl - 1)/(2N C ), N c the number of colours) 



/ 



A 







6C F 






6C F 
N c 





3 

44 
3 

9Vg+4 
3N C 
3(jV c 2 -4) 
N c 



+ 



3 
4 

3N C 



12 


3(iV r 2 -4) 
N c 

-3N r 



(16) 



The non-zero entries of the 4x4 sub-matrices B, C are: Bu = B 2 2 = 8(2N c Nf + l)/3, 
B12 = B 2 \ = 8/3, B 3 i = B 32 = B 4 i = B 42 = 8C F /3, C n = C22 = —26, C 33 = C 44 = —46, 
C 13 = C 24 = 8N c N f /3. The mixing of C 5 into itself is given by 755 = -8(iV c - N f )/3 
Ti|. Note that O5 could mix into O3 at order a s . We find that due to a cancellation of 



different diagrams the corresponding entry 753 vanishes. As a consequence, decouples 
from the mixing at leading order, in agreement with the first reference of |T^ . 

To solve @ with 7(0;^) and P(a s ) evaluated at leading order let 26A i; i — 1 . . .4, be 
the eigenvalues of A and A5 = 7 55 /(26). Let U be the matrix that diagonalises A. Since 
the integration constants must be considered as non-perturbative, we need not keep track 
of factors multiplying these constants, unless they are exactly zero. Thus we only note 
that no element of U vanishes for values of Nf of interest. Since C, and therefore 7^, 
is triangular, we readily obtain 



Ci{a s ) 



4 

E 

k=l 



c 



O' 



1...4 



3 The renormalisation of dimension-six operators has been studied previously in various contexts, 
mainl y Q CD sum rules. We recalculated the entries of 7W in our basis, except for 755, which we take 
from Q. 
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N f 


Ai 


A 2 


A 3 


A 4 


A 5 


3 


0.379 


0.126 


-0.332 


-0.753 





4 


0.487 


0.140 


-0.302 


-0.791 


4/25 


5 


0.630 


0.155 


-0.275 


-0.843 


8/23 


6 


0.817 


0.172 


-0.254 


-0.910 


4/7 



Table 1: Numerical values of A; (N c = 3). 



C 5 (a s ) = Cfa^ 5 (17) 
Q(a s ) = Cf ] a s + £ Cf a7 A * i = 6, 7 



fc=i 
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Ci(a.) = C?a s + cf a? + 2 c£ ] «; Afe < = 8, 9 



k=l 



with a s -independent non-vanishing constants that depend on nine integration con- 
stants and the elements of 7^. The exponents A& are reported in Table |l|. The values 
for Ai to A4 at Nf = 3 are in agreement with [15]. We emphasise again that the coeffi- 
cient functions Cj are independent of the particular Green function or observable under 
consideration. 

We now consider R = q 2 dU(q 2 ) / dq 2 , the 'Adler functions' derived from the correlation 
functions of two vector or axial- vector currents and write their perturbative expansions as 
in ([[]). Having found the a s -dependence of the coefficient functions, we further require 
Roi( a s,q) at leading order. Up to constants, we have RoX a s,q) oc i = 1...4, 
Ro 5 (a s , q) oc a s , Robots, <?) oc aj 1 , i = 6, 7 and RoX a s, q) oc a~ 2 for i = 8,9. Combining 
all factors in (|3|) and taking into account the substitution rule @ to obtain the large- 
order behaviour of R from the imaginary part of I[R] our final result reads 



Ki n 2+x > + K 5 n~ 1+X5 +K 6 + K 8 n 



;i + 0(l/n)) (18) 



for two vector currents.^ The result for two axial currents is identical, except for different 
over- all constants K^. Note that the leading behaviour is generated by the largest eigen- 
value of the anomalous dimension matrix of four-fermion operators. The contribution 
from the three-gluon operator (proportional to K 5 ) is suppressed. The dominant con- 
tribution from this operator comes in fact from its mixing into 06/7 a t next-to-leading 
order (not computed here) rather than the 'direct' contribution through Ro 5 (a s ) in (p~8|) 
above. In general, such next-to- leading order contributions result in l/n-corrections to 
the asymptotic behaviour (P"8"|). Eq. (|T%|) holds when the series is expressed in terms 



4 We display in (|lg) the general structure of the result, although it is not consistent to keep all terms 
without computing 1 /n-corrections to the leading contribution proportional to K\ . 
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of the MS renormalised coupling a s , the convention we assume throughout this note. 
If a different coupling is employed that is related to the MS coupling by a factorially 
divergent series, the coefficients r n change accordingly and (|18|) may not be valid. 

3. The previous result for flavour- U(iV/) singlet currents applies practically unaltered 
to all observables of interest. We briefly discuss them case by case. 



e + e~ annihilation into hadrons. Consider first annihilation through virtual pho- 
tons. The external current is = ip^^Qip, where Qij = diag(e u , e^, . . .) is a matrix in 
flavour space and flavour indices are summed over. Since flavour symmetry is broken 
only by the external current (all quarks are still considered as massless), the 'QCD oper- 
ators' Ci_5 remain unaltered. The basis of 'current operators' Oq-q has to be altered to 
include the operators (trQ) ip'-y^ip d u F Ufl and rf^Qip d v F u ^ instead of Oq. (Similar mod- 
ifications would apply for the axial- vector current.) This ensures that mixing of four- 
fermion operators into the current operators contributes proportional to trQ 2 = J2f e } 
and (trQ) 2 = (J2f e f) 2 , as required by the existence of 'flavour non-singlet' and 'light- 
by-light scattering' terms. The matrices B and C in (|T5|) change, but their pattern of 
non-zero entries does not. Thus, as we are not interested in over-all constants, (fTS]) car- 
ries over to the present case. The leading asymptotic behaviour of the Adler function, 
expressed as in (|TJ), is 

d n n =°° K d flj n\ n 2 ^ 1 ^ = Kd p n n , n i. 97; (1Q) 

where in the last line we have taken Nf = 5. This is to be compared with the large- Nf 
limit [[TIJ which is recovered by setting n 1 - 97 — > n. Eq. (|19|) holds separately for the 



'flavour non-singlet' and 'light-by-light scattering' contributions. In large orders, the 
expansion of the e + e~ hadronic cross section is related to the Adler function by 

e+e~ n^oo d n . . 

r n oc — . (20) 
n 

The suppression by one power of n follows from the fact that the coefficient of the leading 
asymptotic behaviour of the d n is polynomial in the external momentum and therefore 
does not contribute to the discontinuity. 

If we now consider the hadronic width of the Z°, the inclusion of axial currents 
proceeds along the same lines as above for the vector current and modifies over-all 
constants, but again the general pattern of mixing remains the same. The expansion 
coefficients of the Z° width therefore grow as r^ +e in large orders. Note that the fact 
that flavour-singlet terms arise at order a 2 s for axial currents, but for vector currents 
does not affect large-order estimates. 



Hadronic r decay. For flavour non-singlet currents = J £2f= ds V u fU'y fl (l — 75)^/ 
the 'light-by-light scattering' diagrams are absent. For an appropriate basis of current 
operators this is reflected in a change of entries in the matrices B and C. However, their 
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pattern of non-zero entries is not changed. The contour integral that relates the r width 
to the Adler function of the currents suppresses the large-order behaviour just as in case 
of the e + e~ cross section above. We therefore have 

rr«°°< +e -™^n!n a59 . (21) 

In the present case we set Nf — 3. 



Moments of deep-inelastic scattering (DIS) structure functions. The operator 
product expansion allows us to write (at leading twist) 

fidx^FfaQ) =E c *,w(^»^J 4r(A*)- ( 22 ) 

Here the structure function F k can be 2i*\, F 2 /x or F 3 , A l N (ii) denotes (reduced) proton 
matrix elements of twist-2 operators and N (a s ) the coefficients functions. To com- 
pute the large-order behaviour of the perturbative expansion of coefficient functions, we 
need to compute the insertion of dimension-six operators into the quark matrix elements 
of current- current correlation functions.0 Because we consider quark matrix elements 
(rather than vacuum matrix elements as above), insertions of current-current operators 
such as Os/9 vanish and the leading asymptotic behaviour arises from insertion of opera- 
tors such as 06/7- Furthermore, collinear singularities have to be factorised by computing 
the analogues of A l N (fi) between quark states and dividing them out. We choose the MS 
factorisation scheme, in which case the A l N (fi) are pure poles. Dividing them out then 
does not modify the large-order behaviour of the finite terms. Finally, taking moments 
results in a moment-dependent over-all constant, but the n-dependence is the same for 
all moments. Thus, 

CU«.) n =™E K iNWrin 1+M+Xl < +1 . (23) 

n 

Compared to the Adler function in (|19|) one has one power of n less, because the op- 
erators Og/g do not contribute. This reflects that at a given order in a s the diagrams 
that contribute to the coefficient functions have one loop less compared to the diagrams 
that contribute to the Adler function. Eq. (|23|) applies in particular to the case of QCD 
corrections to the GLS and Bjorken sum rule. 



In general, the contribution from IR renormalons to the asymptotic behaviour can 
compete with the UV renormalon asymptotics computed in this note. For the moments of 
DIS structure functions, one finds r^ R /r^ v ~ (— l) n and the nature of the IR renormalon 
singularity is determined by the anomalous dimensions of twist-4 operators [17) . For the 
special case of the GLS sum rule, using the anomalous dimension calculated in [13, we 



J Gluon matrix elements are suppressed and do not contribute to the leading large-order behaviour. 
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obtain 



C GL sM n =°° £/3 "n! [K^ s n ^+M+^ + K G R LS ^i)n n -NPl-{mW c -i/N^ a n+i_ 

n 

(24) 

For the interesting case of Nf > 2, the UV renormalon behaviour formally dominates 
at very large n. For e + e~ annihilation, r and Z° decay the first IR renormalon singu- 
larity occurs at t — — 2/{3 |0| and thus its contribution to the asymptotic behaviour is 
suppressed as r^ R /r^ v ~ (— 2)~ n . Thus, UV renormalons determine the asymptotic be- 
haviour in all cases considered here. However, IR renormalons tend to have large over-all 



normalisation factors in the MS scheme as compared to UV renormalons |L2], Thus, 
dominance of IR renormalons at intermediate n could be expected and is indeed observed 
in the fixed-sign behaviour of the known exact coefficients in the MS scheme. 

The knowledge of the nature of the UV singularity, which we provide in this note, 
could be used to optimise conformal mapping techniques to dispose of UV renormalon 
growth (for an application of conformal mappings in this context, see |^T[ |22[) or to 
improve Pade-type approximations for the Borel transform of the perturbative series by 
combining the information on singularities with the known low-order coefficients. The 
utility of this procedure may be limited since the leading UV asymptotic behaviour con- 
sidered above is not relevant for orders as low as n = 2, the highest order known exactly. 
This is especially so, because the leading behaviour is related to four-fermion operators. 
Subgraphs that contribute to the coefficients of these operators appear first at order a 2 s . 
It takes several more orders to see the exponentiated effect of two-loop running of the 
coupling and the eigenvalues of the anomalous dimension matrix embodied in the factor 

n /3i//32+A!_ 
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